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We compute the free energy density of QCD to order as at very high density and 
^ \ non-zero quark masses. The counterterms needed to renormalise the theory to order is 
same as the vacuum (non-zero density) theory. We investigate the response of the theory 
^ . to non-zero quark masses and chemical potentials. We study quark number density and 
■ quark number susceptibility in the high density limit, where the ratio of the quark mass to 
^ ' the corresponding chemical potential is very small (m/// <^ 1). In this limit both number 
^ ' density and susceptibility contain a ln(m//i) contribution at order ag- We compute the 
scalar and pseudoscalar susceptibilities to order in three flavour QCD at high density 
'-^ \ taking quark masses and chemical potentials to be degenerate and non-zero. At extremely 
^ . high density, since ag is very small, both the susceptibilities are found to be same in the 
^ ■ chiral limit. This means that the scalar-pseudoscalar splitting is absent in the CFL phase. 

> 

S ■ 1 Introduction 

The study of QCD at very high densities and low temperatures has been received a lot 
of interest in the last few years. It has been argued that at a sufficiently high density the 
colour symmetry is spontaneously broken, leading to colour superconductivity 0, ^ ^ 
(see for recent review). This novel phase might be accessible in the core of neutron 
stars or in the baryon-rich fragmentation region of heavy-ion collisions. 

At extremely high density, where the density is of the order of 3 to 5 times the 
saturation density of nuclear matter, the wave function of the quarks in nucleons will 
overlap with that of quarks in other nucleons due to asymptotic freedom. At such high 
density quarks are no longer confined in nucleons and thus the nuclear matter will become 
a quark matter . The matter at such high density will consist of a Fermi sea of essentially 
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free quarks and its behaviour is dominated by the high momentum quarks that hve at the 
Fermi surface. It was shown by Bardeen Cooper and Schriffer (BCS) |^ that in presence of 
attractive interactions a Fermi surface is unstable. Since in QCD the attractive interaction 
is provided by one-gluon exchange between quarks in a colour antisymmetric channel, the 
true ground state of the system is not the naked Fermi surface, but rather a complicated 
coherent state of particle and hole pairs or condensates, called "Cooper pairs". The 
formation of condensates leads to energy gaps for both the (quasi) quarks and some 
or all gluons. At asymptotically high density, and because of asymptotic freedom of 
QCD, it is possible to compute the quark gap and the gluon masses from first principles 
i, I, 

In order to obtain the bulk quantities of a system, such as number density of particles, 
susceptibilities etc., it is necessary to obtain the free energy density or effective potential 
of the system. These bulk quantities give us informations about the behaviour of the 
system in the microscopic level. Freedman and McLerran computed the renormalisation 
group (RG) improved thermodynamic potential at finite density to order in non-abelian 
gauge theory with massless quraks [|T^ . Using this RG improved thermodynamic potential 
they obtained the effective quark number density. They also obtained using RG equation 
(Gell-Mann-Low equation) the flow of strong coupling constant as a function of chemical 
potential. Their evaluation shows that as decreases with the increase in density. The two- 
loop exchange contribution in the thermodynamic potential for massive quarks, later, has 



been evaluated in the ref.[|T5 



In the present work we have studied the susceptibilities to leading order in strong 
coupling constant (order a^) in the high density phase of QCD. As a first step of this we 
have computed the free energy density to order at finite density taking quark masses 
to be non-zero. For massive quarks the cancellation of ultraviolet (UV) divergences is 
a bit involved issue. The free energy at order is UV divergent. This UV divergent 
piece has two parts : One is the vacuum (zero density) part and the other one is non-zero 
density part. However, both the parts get cancelled from the counterterm diagrams. The 
counterterms needed to make the free energy density finite are the same as what we had 
in the vacuum theory. Using this renormalized effective action we have computed the 
quark number density and the quark number susceptibility to order a^. We find that 
at finite baryon density, which is high enough too, the density of states of a particular 
quark flavour near the Fermi surface decreases with the increase in mass of the quark. 
Therefore, at finite baryon chemical potential the density of states of strange quark is 
much lesser than that of u and d quarks near the Fermi surface. 

The colour and flavour symmetry of QCD with three massless quark flavours break 
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at high density. The breaking of this flavour and colour symmetry leads to eight massive 
gluons and gapped (quasi) quarks. The nine (quasi) quarks (three colours times three 
flavours) fall into an 8 © 1 of unbroken global SU(3). Therefore, in order to parametrise 
this phase completely one needs two gap parameters. This particular phase is called 
the colour-flavour locked (CFL) phase [Q. In the CFL phase both colour and flavour 
symmetries are dynamically broken by the formation of quark-quark condensate near the 
Fermi surface. Since the chiral symmetry in this phase is broken as well by the formation 
of condensates between quark and antiquark pair, the phase can be best described by an 
effective theory of pseudo-Goldstone bosons 0, |T^. Since the scalar and pseudoscalar 
susceptibilities are related to the order parameter of the chiral symmetry breaking {{i/jiIj), 
where ip represents the quark field) , it is pertinent to study both the susceptibilities in 
the CFL phase. Using the renormalized free energy density we have computed the scalar 
and pseudoscalar susceptibilities to order in the high density phase of QCD with three 
flavours, taking quark masses and chemical potentials to be degenerate and non-zero. We 
find that at extremely high density, since is very small, both the susceptibilities are 
same. 

We organize the paper as follows. In sec. 2 we compute the renormalized free energy 
density to order a^. In sec. 3 we compute the susceptibilities. In sec. 4 we discuss the 
results. In Appendix-I we give the fermion propagator in the dense medium. In Appendix- 
II, III and IV we discuss the estimation of integrals which are needed to compute the free 
energy density in sec. 2. 



2 Computation of free energy density to order as 

The gauge Fixed QCD Lagrangian with Uf number of quark flavours at finite density is 
given as 

+9^ci(rt^ - ^r^M^'^)Q„ (2.1) 

where F^^, = — S^AJ^ — gf°-^'^A'^^A^^. a,b,c (= 1, ■ ■ ■ 8) are the gauge group indices 

and A,B,C {= 1,- ■ - nj) are the flavour indices of quarks. ^ is the gauge fixing parameter 
and we shall take ^ = 1 (Feynman gauge) in the following calculations. 
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2.1 Counterterms to order 

At order the QCD Lagrangian is divergent. The divergences are coming from the one- 
loop momentum integrals in the quark and gluon two point functions. In order to remove 
the divergences one needs counterterms. The renormalized QCD Lagrangian reads 

In order to obtain the counterterm Lagrangian let us first calculate the quark two point 
function at order g"^ (the quark self energy diagram is shown in Fig. la). We calculate the 
self energy diagram in Fig. la using the quark propagator in Appendix-I. This diagram is 
the sum of two parts: The vacuum part and the finite density part. Since the vacuum 
part is ultraviolet (UV) finite, we shall calculate only the vacuum part of the quark self 
energy. The vacuum part of the self energy reads 

S(^"'^)(p)ab = -g^C2{R)5AB^^' J ^A^^iSip - k)Y]iD,.{k), (2.3) 

where the gluon propagator iD^,^{k) — — (in Feynman gauge), A is the regularisation 
scale and d (= 4 — 2e, e > 0) is the space-time dimension. Here C2{R) — | (for SU(3)) 
and p — {po + iiA,p)- Performing the k integration in d dimension we obtain 

r(e) f',,, - e)tp - {2 - e)mA 



^'-\p)ab{p) = ^9'5AB{-ir^^^^^^, £ dt{l - t)-H 



67r2(47rA2)-e Jo ^ ' _ :<^y 

(2.4) 

This term is divergent when e — > 0. In order to remove this divergence we add to C the 
following counterterm (the corresponding counterterm diagram is shown in Fig. lb) 

- - Q^mL^yA -^rm / dtr^ii - tri^AOi., (2.5) 



where 



= 



67r2(47rA2)- 
1 - e)t{^^{ido + ^^A) - - m^(2 - e) 



(2.6) 



For minimal subtraction we shall retain terms up to zeroth order in e. 

Next we calculate the gluon two point function at order g"^ (the gluon polarisation 
diagram is shown in Fig. 2 a). Although the diagrams with gluons and the ghosts inside 
the loops give a non-zero contribution to the gluon polarisation tensor at order g^^, they 
contribute zero to the free energy. So for our purpose we calculate the gluon polarisation 
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taking only quarks in the loop. The vacuum part of the gluon polarisation (since the 
finite density part is UV finite) reads 



(27r)^ 



trs[iS{k)YiS{k-p)Yl 



(2.7) 



where iVg means the trace of the product of Dirac gamma-matrices. We evaluate the 
traces in d space-time dimension using the formulae 7^7^ = d and trsij^ju) — f{d)T]^v 
{f{d) ^ 4 as 0? ^ 4). We drop the terms linear in k and obtain 



Jo 



dtk^^- — -T 

{2'nY 



2k^k^ - rji^^k^ + r]^,^{m\ + pH{l - t)} - 2p^p^t{l - t) 

[P - {m\-pH{l-t)}]' ' ^ ' ' 

where wc have used the Feynman parametrisation. Performing the /c-integration in d 
dimension we obtain 

r(e) 



167r2(47rA2)- 



C dtrHi-t)- 

Jo 



X 



t{i-tx 



t(\-t) 



-1+e 



-I- 



m\'>l^l<^ -PnPut{l - t) 



2 

P t(i-t) 



(2.9) 



This term is divergent as e 0. In order to remove this divergence we add the following 
counter term to the Lagrangian (the corresponding counterterm diagram is shown in 
Fig.2b). 

C?? = - /(d)(-l)-T(e) dtr^l - t)-'AlO^"'Al, (2.10) 

Jo ^ 

where 

- t) _ m\if " + f (1 - f^di^d" 



327r2(47rA2)- 



-□ - 



t{l-t) 



-1+e 



+ 



-□ 



t(l-t) 



(2.11) 



For minimal subtraction we shall retain terms in upto zeroth order in e. Therefore 
the renormalized Lagrangian to order g'^ is 



(2.12) 



where C^^t is a field independent counterterm which is required to renormalise the free 
(gf-independent) part of the free energy density. This can be written as 

£2^-C(m^,e), (2.13) 

where C is a functions of quark masses (m^) and e. For minimal subtraction we shall 
retain in C the terms upto first order in -. 



5 



2.2 Free energy density to order 

The renormalized QCD partition function reads 

^ren = y 2^0e^^'^'"^'^\ (2-14) 

where the variable set 4> = {'4^, V', c^, Ca}- The renormahzed free energy density is 

F=^lnZren: (2.15) 
III 

where Q is the 4-space-time volume. Retaining only the gi-independent quadratic terms 
in Sren wc expand the remainder up to order in perturbation theory and integrate 
over 0. The graphs that will contribute to order g'^ are shown in Fig.3. Therefore the 
renormalized partition function to order g'^ reads 



^ren ^ren 



(2.16) 



where /^^''^ is the contribution of the diagram in Fig.3j (j = a, b, c, d, e, f). Here Zf^i reads 

Uf 

Zll = e^"' n {det[i^ - + ix^l' - e^'d^] }" , (2.17) 



where the power 3 in the determinant is for three colours associated with each flavour 
of quarks. The evaluation of this determinant is discussed in the Appendix-II. The final 
result is 



exp 



3m^ rs , ( m\\ 



+ 



4^2 ^--v-A -Ayr-A 2-"^y ' 2 ^ V mA 

where we have chosen C = J2^=i ^'^^ minimal subtraction. 

We write in the following the contribution of the diagrams in Fig.3. 



,(2.18) 



Ia 



(27r)4 (27r)4 (p - q)^ 
(3b) _ 2 I d^p dS A^'^'A^,, 



^ J (27r)4(27r)4pV(p + g)'' 
^ ^ i (27r)4(27r)>V(p-g)'' 
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r{3d) 



r(3e) 



r(3f) 



d^p d^q B 



(27r)4(27r)4pV' 

i9\-l)- ' 
27r2(47rA2)-- jo 



r(e) j^dtr%i-tr^ j ^,trsp{k)iS{k)l 



(27r) 



47r2(47rA2) 



:^J{d)m f^dtt-^l-trj -^,[0^''{k)iD,^{k)l (2.19) 



where 



B 
0{k) 

0'"'{k) 



= ri'^'ip - qY + r/'^"(p + 2qY - r/'^^(2p + g)^ 
(1 - e)t^-mAi2 - e) 



fc2 



77'^'^t(l - t) m\ri'"' - t{l - t)k''k'' 



k^ - 



t{i-t) 



-l+e 



+ 



t(l-t) 



(2.20) 



Only the diagrams in Fig. 3a, 3e and 3f will give non-zero contribution to the free energy, 



the rest will not contribute to the free energy i.e. 



(3b) 



r(3c) 



r{3d) 



0. The integral 



r(3a) 



can be written as the sum of three terms : The vacuum part, density dependent UV 



divergent part and density dependent finite part. So we write it as 

/r = V[/(^)/r + /(^)/? + /f], 



(2.21) 



where we have evaluated the trace over the Dirac gamma matrices in d space-time dimen- 
sion and write the d-dimensional integrals in the following: 



r(0) 



r(i) 



A 



4e 



d'^k d'^p 



(2 - d)k.p + m\d 



7-(2) _ 



{2'kY (27r)'^ (F - mi + ie){p^ -m\ + ie){{p - kf + ie} ' 

(2 - d)k.p + m\d 
(/c^ — 772^ + ie){p'^ — m\ -\- ie){{p — kY + ie} ' 

d'^k d'^p 



X 



-167r 



X 



2m\ 



{2ttY (27r)4 
k.p 



5{k' - m\)5{p' - m\)e(ko)e{po)9{iiA - ko)e(iiA - Po) 



m. 



The integral i^^®-* can be splitted up into two parts: The vacuum part and the density 
dependent UV divergent part. Going over to d space-time dimension and evaluating the 



k.p 



(2.22) 
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spinor traces we can write it as 



(2.23) 



where 



t(0) 



t(1) 



^(-1)- 

(47rA2)- 



-r(e) £ dtt-'{l-t)-'A''' j 



dfik (1 - e)tk'^ - m\{2 - e) 



{k"^ — m\ + ie) 



1^2 '"A 



UttA^ ) 1-26 



2 - 2e 



-(2-6) 



X 



27r(5(r - m\)9{ko)e{iiA - ko) 



In a similar fashion we can write going over to d space-time dimension as 



(2.24) 



7-(3f) _ 

-'A — 



47r2 



where 



(2.25) 



d'^k 



(27r)' 



t(i-t) 



(2.26) 



In the minimal subtraction scheme we shall retain terms in J^2\ J^^ and K^^ upto first 
order in K The evaluation of the integrals J^^ and J^^ are given in the 

Appendix (III and IV). The results are written in the following. 



r(i) 



fm^\ (my 



-2e 



3 7 67s T . , T . « 2 

— J - - + 15 + lA-iE - 67e - — 

e e 2 



6477^ 



■niA 



7-(2) _ 



1 



3m, 



27r4 



In 



+ 



IGtt^ 



In^ 



ruA 



m\{nA - mAY_ / V/iA + + VA«a - mA 



47r4 



2A/iUA 



8 



(0) ^ rn^ f m\ \ 



-2e 



J 



7(1) _ 

'J A — 



K 



(0) 



327r2 I 


,47rA2 j 








V47rA2 


mi / 




327r2 \ 


, 47rA2 j 



[3 7 6 1 

+ 7e 



fJ'A 



rriA 



-2e 



3 7 6 
-J + - - -7£ 



(2.27) 



The vacuum UV divergences in the diagram Fig. 3a get cancelled from the counterterm 
diagrams Fig.3e and 3f and the density dependent UV divergence in Fig. 3a gets cancelled 
from the counterterm diagram Fig.3e. Therefore the renormalized free energy density to 
order as reads 

F^^fA, (2.28) 



A=l 



where 



3m\ 



3 , { m\\ 



+ ^^^^A/fA-rn^ (/^A - l^A 



-m\ In 



rriA 



+- 



7J-3 



(15 - 147^ + 671 + y) - kl^l - Oif^l - 20 



(2 - ^7e)"^>aV/4-^+ 



rriA 



+m^{ (2 - --fE)m\ + ^^^A\/ fiA - m\ - 4/x^} In 



^^A + V/^A - ^1 



rriA 



-2m\{iiA - mA)^ln 



V/M + mI+ V/iA - niA 



+ 0(a 



(2.29) 



2^/^A 

3 Response to chemical potentials and non-zero quark 

masses 

In this section we shall investigate the response of the system with the change in chemical 
potentials and quark masses. The quark number density is obtained by taking a derivative 
of the free energy density with respect to chemical potential. In the high density limit 
this is given by 

dF 



Ua = 



9 



2 fx A, 



+0(a„— )• (3.1) 

Here we have retained terms upto zero-th order in fiA- In the chiral hmit with non- 
interacting quarks the chemical potential /xa ~ ^A- However in the non-chiral limit with 
interacting quarks the chemical potential cannot take such a simple form in terms of 
density. 

The quark number susceptibility to order in the high density limit (^ <^ 1) reads 



XA 



dfiA 

~~^I^A ~ 7r~^^A "I ^ 



-6fil + 3(4 + ^eWa - 14mi In 



+0{al—), (3.2) 

fJ'A 

where we have retained terms upto zero-th order in /i^. The susceptibility has a simple 
physical interpretation. It measures the density of states of a particular quark flavour 
near the Fermi surface. At finite baryon density {fiA = for all A, where /2 is the baryon 
chemical potential) as the mass of the fermion increases the susceptibility or the density 
of states near the Fermi surface decreases. Therefore the density of states of the strange 
quark near the Fermi surface at finite baryon density is smaller compared to the u and d 
quarks. 

The colour and flavour symmetry of QCD with three massless quark flavours break 
at high density. This breaking leads to eight massive gluons and gapped (quasi) quarks. 
This particular phase is called CFL phase ||^. We shall calculate in this phase the scalar 
and pseudoscalar susceptibilities which are related to the chiral order parameter of the 
theory. In order to do this we consider the three flavour QCD with u, d and s quarks 
at high density. We take quark masses and chemical potentials to be degenerate; m„ = 
rrid = rris = m and Hu = f^d = f^s = /W- This means that the theory respect the SU (3) l+r 
flavour symmetry. Let us define the following scalar and pseudoscalar correlations: 



Cs{q) = J rf^xe^^-^'(0 I ^{x)iIj{x)^{0)iIj{0) \ 0), 

Cp{q) = J d^a;e^''-^(0 | ^{x)i-f5^{x)^{0)i-r5^{0) \ 0). (3.3) 
The scalar and pseudoscalar susceptibilities are defined as, 

Xs = Cs{0) and xp = Cp{0). (3.4) 
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These two susceptibilities are related to the chiral order parameter in the following way 

d{il>il)) 



Xs 



and Xp 



(3.5) 



dm ~ m 

where = {uu + dd + ss). The chiral condensate is obtained by differentiating the 

free energy density with respect to m. Once we obtain (ipip) we can compute the scalar 
and pseudoscalar susceptibilities using eqn.(^.5|). Both the scalar and pseudoscalar sus- 
ceptibilities in the high density limit read 



Xs 



9m 



2 r 



27r2 
3a, 



+ 



— S'jE — 3 In 

;i3 - 37£;)yU^ 



f m^ \ 
1^47rA2 ) 

19 - 



3 



2 o 2 9 2 1 / ^ 

u — 6m m m — 

^ 2 \2^^^ 



39 

y 



lE + 971; + 7^ \m 



4 



m 



m 



m 



+ (187s - h)m' In — + In^ — + In — 



9m 



2 r 



Xp 



27r2 
3a., 



7£; 



In 



2/i 
I m^ \ 



+ 



2 



m In 



m 



+ 



(8 - 37£;)/x' - |lO - y7i^ + ^ll + ^vr^j 



m 



67£;)m^ In 



m 



+ 2m2 In^ 



m 



+ 2/1^ In 



m 



(3.6) 



At very high density (^ <^ 1) as we take the chiral limit (m 0) both the susceptibilities 
are logarithmically divergent at order a^, i.e. xs ~ — ^{^^ ~ ■^Ie) + ^ In (^) | 

and Xp ~ + v(8-37i?) + ^ln (|^)|. Since at very high density a, < 1 0, the 

resulting contribution at order as is small compared to the leading order (zero-th order in 
as) one. Therefore, at extremely high density both scalar and pseudoscalar susceptibilities 
are same in the chiral limit, i.e. Xs = Xp = ^Z^^- 



4 Discussion and Conclusions 

We have computed the free energy density of QCD to order a, at finite density taking 
quark masses to be non-zero. The free energy density is found to UV divergent. The 
counterterms needed to renormalise it is the same as what we had in the vacuum theory. 
We have computed the quark number density and the quark number susceptibility in the 
high density limit {mA/fJ'A 1) using this free energy density. In this limit both of them 
contain a ln(myi/yU^) contribution at order a^. 
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We have computed the scalar and pseudoscalar susceptibihty in the high density phase 
of three flavour QCD taking masses and chemical potentials of the quarks to be degenerate 
and non-zero. We find that at very high density, since -C 1, both the susceptibilities 
are same in the chiral limit. Since the mass is related to the susceptibility via = Zx~^ 



T9[] , where Z is the wavefunction renormalisation, and since the scalars and pseudoscalars 
renormalise in the same way, we obtain the following equality in the chiral limit: M| = 
Mp. Therefore in the three flavour QCD both the scalar and pseudoscalar masses are 
degenerate in the CFL phase. 

Acknowledgements: I would like to thank Debadesh Bandopadhyay, Subhasis Basak, 
Amit Kundu and Pushan Majumdar for useful discussions. 



Appendix-I: Fermion propagator at finite density 



Free fermion propagator at finite density reads 



{2tt)^ ^ — m — yU7° + ze7°po 
d'^P -fp.x 7°(Po + /^) -l-P + m 



where Up = \/p^ + ni^ . Performing the po integration we obtain 

f ^l+^nf^_ ,.-ip.x+^^xO 



J (27r)'^ L 2po 2po 
where Pq = ujp. After doing some algebra we can write it as 



where iS{p) = ij) + m)iS{p) and 



(i-i: 



(1.2) 



tSip) = - 2nS{p' - myipo)9ifi - po). (1.4) 



Appendix-II: Evaluation of Z^^2% 
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Define 

D^i^-mA + I^Al° - e7°9o. (II.l) 

Now 

detL> = e^'''"^, (II.2) 
where Tr is taken over spin matrices and space-time. 



TrlnD 



= J d'^x{x I tr,]nD \ x). (II.3) 



Inserting complete set of eigenstatcs of momentum in between and using the formulae 
{k\x) = e-''^-^, {k' I k) = {27ry6\k - k') we obtain 

r d^k 

TrlnD = J d*Xj^trslnD{k). (II.4) 

Here we have used the fact that operator D is diagonal in the momentum representation 
i.e. {k\D\k') = {27i)^6^{k - k')D{k). Using the formula trluD = Indeti:) we obtain 

TrlnD ^ [ d^x^^\ndetsD{k), (II.5) 
J (2'?r) 

where the det^ means the determinant of the resulting spin matrix. Evaluating the deter- 
minant of the spin matrix we obtain 

TrlnD ^2 [ d^x^^ ln{k^ - m\) , (II.6) 

J (2'7r) 

where k — {ko + /ia + iekg, k). Therefore the renormalized free partition function reads, 

[I 

A=l 

e^^' exp jei] i: / ^ ln{F - m^) | . (II.7) 



We use the following parametrisation 



In A;' - ml = / — + / — , II.8 

^ ^' kQ + {l-%e)x Jo ko + {l-ie)x ^ ' 



where cj^^ = \Jk'^ + m\. Here we have ignored the mass(m^) and density(/i^) indepen- 
dent terms. We first do the /co-integration going over to the d space time dimension and 
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then do the integration over x. Next we do the remaining k integrations in c? — 1 and 3 
dimension for the vacuum and the finite density part respectively. The result is 



A 



2e 



167r2 

i 



1 3 
7+2 



lE 



In 



247r2 

3 
2 



( 

1 ^2 
- 2^A 



, - 4 1 / + \/^^A - ml 
+ -m^ m 



(11.9) 



Appendix-III: Evaluation of I^^ and jjf^ 



Performing the momentum integration over k and p the integral /jf^ reads 



r{0) 



2+e 



x|(2?/-l)^ dxx~'{l- xyy^^' + (l-y) dxx~'{l- xyY 



3+e 



Define 



In = C dxx~\l - xyy (n>0), 
where it satisfies the following recursion relation 

1 



1 - e 



Using this recursion relation we obtain, 



'2-e 



-l+e 



Therefore we can write 



r(0) 



1-e 



'3-e 



-2+e 



-l+e 



2-e 



(l-e)(2-e) 



(£^)"^2-e)r(-l + 2e)f... 



-2+e 



V167r2^ 

x{(2|/-l)/2_, + (l-y)/3-e}. 



Using eqn.( [lll.^) and performing the integration over y we obtain 



r(0) 



m. 



f rnA \ 



-2e 



167r2 / V47rA2; 



7 6 9 

- + -7i5 - 15 + 147s - 67I - — 



.(III.l) 
(III.2) 

(III.3) 
(III.4) 



(III.5) 



(III.6) 
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Performing the integration over p the integral I^^ reads 



r(l) 



X 



(27r)32a;fcA 



m\ ( m\ Y 
647r4 V47rA2 J 



4 + 37^ 



X 



^^A + - "^A 



rriA 



where 



(27r)32a;fc^ 



87r2 



I^A\/l^A-m\ - m\ In 



flA+ yj^A^^A 

rriA 



(111.7) 



(III.8) 



The integral /^^^ reads 



r(2) 



1 



HA+ \/jA^2^A 

ruA 



(III.9) 



where 



La — 2//^ 



(1 — Q;)^lnQ;+ dx dy lln \x — y\ —In \xy + Vx^'^^o^^fv'^ 

J a J a 



'x^ — a^yy^ — — a \ 

(III.IO) 

and a — m^/AfA- We can carry out the first integral and the result is 

r dx t dy\ii\x-y\ = (1 - a)^ ln(l - a) - hl-af. (III.ll) 

J a J a 2 

In order to do the second integral we change the variables from x and y to 9 and (f) 
respectively where x — a cosh 9 and y — a cosh 0. The integral becomes 



j dx I dyhi \xy — y, 

J a J a 



x"^ — oP'sjy'^ — oP- — o?\ 



= 2{l-af \a a +{l-af 1^2 + 20^ dicoshO) d(cosh(f)) 

Jo Jo 

xlnsinh l^^^j , (111.12) 

where u = In ^^^^ -^^ ■ The remaining integration over 9 and can be done easily. 
Finally the result for La is 



La = ii\-m\- In 



/^A + ^f^A - ^A 



rriA 
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ruA 



+ m\ In^ 



' llA + ^fJ'A-ITT'A^ 



rriA 



-2{iJ,A - ruAfln 



y/fiA + niX+ s/liA - rriA 
2^ 



(111.13) 



Appendix-IV: Evaluation of Ja^ and Ja^ 



Performing the integration over k the integral Ja^ reads 



= -1^5 (4^)"" r(-2 + e)(l - .)(2 - - A-J, (IV.l) 



where 



Jo Jo [L — 

X, = 2 ['dyy-%l-y)-'+' ['dx— 
Jo Jo (1 - 



X 



-1+e 



y + xi/) ^+^^ 



-l+e 



(IV.2) 
(IV.3) 



y + xy) ^+^^ 

Let us first evaluate Xi. Performing a integration by parts with respect to x we obtain, 



^1 



-l+e 



1 2(1 



e e(l + e)' 



+ 



1 -e)(l -2e) p /■! , 1 
-—l^—^y^ dxx'+'{l-y + xyy 
e(l + ej ./o 



2e 



(IV.4) 



There remains a integration over x within the third bracket. This integrand is finite at 
e = 0. So we first expand the integrand about e = and then carry out the integration 
over X. 



Xi = fdyy-\l-y) 
Jo 



-l+e 



1 2(1 -e) 2(l-e)(l-2e) rl 



.2 I ^,„.2 



e e(l + e) e(l + e) 

6(1 - yfHl -y) + ^e^y^ - \e'y + • • • + 0(6^)} 



(IV.5) 



Here dots imply the terms of order e^, which after integration over y will contribute to 
order e in Xi. Since we are using minimal subtraction scheme, those terms of order e in 
Xi are not relevant for us. Carrying out the integration over y and retaining terms upto 
zeroth order in e we obtain 

^1 = ^ + 1 (IV.6) 
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In a similar fashion we can carry out the integration in X2 and the result is 

X2 = ^ + 4. (IV. 7) 



Using the result of Xi and X2 we can obtain J\ 



(0) 



Using eqn. (|III.8|) we perform the integration over k and the result is 



(1) ^ 3mj^ / ruA y 
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Figure Captions 



Figure 1 : (a) Self energy diagram of quark at order ag- (b) The counterterm diagram 
of quark two-point function at order a^. 

Figure 2: (a) Gluon polarisation diagram at order with only the quarks inside the 
loop, (b) The counterterm diagram of gluon two-point function at order ag. 

Figure 3: Diagrams contributing to the free energy density at order Ug- Here solid 
lines represent the quarks, coiled lines represent the gluons and dashed lines represent the 
ghosts. 
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Figure 2. 



20 




Figure 3. 
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